Diffusion and mixing in gravity-driven dense granular flows 
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We study the transport properties of particles draining from a silo using imaging and direct 
particle tracking. The particle displacements show a universal transition from super-diffusion to 
normal diffusion, as a function of the distance fallen, independent of the flow speed. In the super- 
diffusive (but sub-ballistic) regime, which occurs before a particle falls through its diameter, the 
displacements have fat-tailed and anisotropic distributions. In the diffusive regime, we observe very 
slow cage breaking and Peclet numbers of order 100, contrary to the only previous microscopic model 
(based on diffusing voids). Overall, our experiments show that diffusion and mixing are dominated 
by geometry, consistent with fluctuating contact networks but not thermal collisions, as in normal 
fluids. 
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Granular flow is an attractively simple and yet sur- 
prisingly complex subject Fast, dilute flows are 
known to obey classical hydrodynamics (with inelastic 
collisions), but slow, dense flows pose a considerable chal- 
lenge to theorists, due to many-body interactions and 
non-thermal fluctuations. Beyond their fundamental sci- 
entific interest, such flows have important engineering ap- 
plications 0, e.g. to new pebble-bed nuclear reactors 
whose efficiency and safety depend on the degree of mix- 
ing in very slow granular drainage (< 1 pebble/min). 

Although dense granular drainage is very familiar 
(e.g. sand in an hourglass), it is far from fully under- 
stood. Over the past forty years, a number of theo- 
retical approaches have been proposed for steady state 
flow 0, @j 13 ■ Continuum approaches are based on 
the critical-state theory of soil mechanics and yield only 
mean velocity fields 0il3- On the other hand, the diffu- 
sive void model 0, takes a particle approach, in which 
'voids' injected at the orifice cause drainage by diffus- 
ing upward and exchanging position with particles along 
the way. Averaging over the void trajectories yields the 
same continuum velocity field for particles as the 'kine- 
matic model' 0, 0, which provides a reasonable fit to 
experimental data with only one fitting parameter (the 
diffusion length, b) 5, 6, 9, 10], although the void model 
on a regular lattice (as in Ref. |^) underpredicts its 
value [l^. Remarkably, these models depend only on 
geometry and not on momentum, energy, etc. 

In spite of the success of the kinematic model, how- 
ever, its only microscopic basis, the void model, greatly 
overpredicts diffusion. To see this, consider the Peclet 
number, PCj. — YyCL/D^, the dimensionless ratio of ad- 
vection in uniform downward flow of speed, v^ , to diffu- 
sion with a horizontal diffusivity, D^, at the scale of a 
particle diameter, d. In the void model, when a parti- 
cle falls a distance Aft,, Ve^ = {Ah/At)d/{{Ax'^)/2At) = 
Ah2d/{Ax'^), is of order one for any conceivable packing 
since Ah « Ax « d, and therefore it diffuses horizontally 
by roughly V Ah. This prediction is contradicted by ev- 




FIG. 1: An initially flat, off-center interface between two 
regions of differently colored beads (a) stretches and roughens 
after draining half of the silo (b), but little mixing is observed, 
(c) Contour plot of the average downward velocity fleld, Vz, 
with an orifice width, W = 16 mm = 5.3 d, where d is the 
particle diameter. The white box indicates a region of nearly 
uniform flow where all subsequent measurements are made. 



eryday experience and our experiments below, which ex- 
hibit far less mixing. An attempt to resolve this paradox 
with a new model appears in a companion paper |l2l | . 

In this Letter, we describe particle-tracking experi- 
ments on silo drainage using similar techniques as in a 
recent (lower-resolution) study of the velocity field [isj . 
We focus on the statistical evolution of particle displace- 
ments and topological "cages", which should aid in de- 
veloping new microscopic models. Our data may also 
have implications for recent attempts to apply thermo- 
dynamic approaches from glassy dynamics to granular 
flows [13) fl3 flfij l . Although we do not deflne a "granular 
temperature" , we observe the presumably related effect 
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of varying the flow rate in all of our measurements. 

Our experimental apparatus involves glass beads {d = 
3.0 ±0.1 mm) in a quasi- two dimensional silo (20.0 x 
90.0 X 2.5 cm = 67 X 300 x 8.3 d). The particles are 
observed near the front wall made of transparent glass, 
where the slight polydispersity reduces the tendency for 
hexagonal packing. (As seen in Fig. ^ there is no long- 
range order, although the wall induces some short-range 
order that may affect our results.) We find that varying 
the thickness of the silo has insignificant effect on the 
diffusion properties and therefore we report our data for 
a single thickness. We track individual particles using 
a high-speed digital camera with a maximum resolution 
of 512 X 1280 pixels at 1000 frames/sec. Particle posi- 
tions are obtained to sub-pixel accuracy using a centroid 
technique (±0.003(i, d ~ 15 pixels). 

Our first experiment provides a visual demonstration 
that particles mix much less than predicted by the void 
model. We load the silo with black and white (but 
otherwise identical) glass beads, forming two separate 
columns, as shown Fig.^a). From Fig. ^b), where half 
of the particles have drained (in 30 sec), it is clear that 
the black-white interface has not smeared significantly, 
although it has roughened. The small degree of mixing 
is consistent with the segregation of bi-disperse beads in 
a similar apparatus [lol |. 

The mean downward velocity field, Vz, is shown in 
Fig. n^c) throughout the silo, as measured by direct par- 
ticle trackin g. A distributed filling procedure similar to 
that in Ref. [llf was used to add grains to the silo. Then 
the orifice was opened and steady state flow was allowed 
to develop before acquiring the data used for subsequent 
analysis. The flow speed is highest at the center-line, 
with a maximum near the orifice, and decays to zero to- 
ward the sides. Fairly good agreement with the kinematic 
model is obtained with b — 1.3d. 

In order to investigate particle dynamics in a simple 
setting, we focus on a small 17 x 87 d region of nearly 
uniform flow, the white box in Fig. ^c) . For the mea- 
surements below, we track about 1370 particles through 
this window for 4 seconds in steps of 1 ms. The average 
flow speed v^ , the only control parameter in this study, is 
varied by changing the width of the orifice, W. The flow 
is fairly smooth for > 8 mm (about 3d) and nearly 
continuous for W > 16 mm. For simplicity, we vary W 
in the range 8 mm < < 32 mm, in increments of 4 mm, 
to avoid the complicated regime of intermittent flow . 
This corresponds to 1.38 d/sec < v^ < 18.39 d/sec, and 
data is consistent with v^ cx {W — d)^'^. We compile 
statistics from all tracked particles in six experiments per 
flow speed (except two for W = 32mm). 

From the positions of the particles sampled at 1 ms in- 
tervals, we calculate the horizontal and vertical displace- 
ments. Ace and Az, relative to a frame moving with the 
mean speed of the flow. A typical trajectory computed 
in this way in Fig. Ola) shows periods of small fluctua- 
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FIG. 2: (a) A typical particle trajectory sampled at 1 ms 
intervals in a frame moving with the average flow speed, v^. 
(b)-(c) Normalized PDFs for the 1 ms particle displacements, 
Ax and Az, for various flow speeds, v^, compared to a stan- 
dard Gaussian distribution (dotted line); standard deviations, 
ax and a^, are of order 10~^d. (d) The kurtosis of Aa; versus 
the mean distance fallen, v^At. 



tions with occasional, much larger steps. This suggests 
that the probability density functions (PDFs) for Ax and 
Az (for At — 1 ms) should have fat tails compared to a 
Gaussian, which is confirmed in Fig. |2Ib)-(c). 

Fat-tailed PDFs have also been observed in colloidal 
glasses and attributed to cage-breaking but a spe- 
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FIG. 3: (a)-(b) Mean squared horizontal ((Aa;^)) and vertical 
((Az^)) displacements versus mean distance dropped, which 
collapse onto a single curve for different flow speeds, v^, ex- 
cept for the smallest where the flow is intermittent, (c)-(d) 
Diffusion coefficients (D) and Peclet numbers (Pe ) in the hor- 
izontal (x) and vertical {z) directions versus Vz. 



cial feature here is the asymmetry of the PDF for Az in 
Fig. |2tc). Downward fluctuations (Az < 0) are larger 
than both upward (Az > 0) and horizontal (Aa;) fluc- 
tuations. We attribute this to the fact that particles 
are accelerated downward by gravity while being scat- 
tered in other directions by dissipative interactions with 
neighbors. 

Looking again at Fig.|2Ia), it seems that the large fluc- 
tuations in particle displacements would be reduced by 
coarse-graining in time, perhaps enough to recover stan- 
dard Gaussian statistics. Indeed, as shown in Fig.[2Id), 
the normalized kurtosis, k.^ = (A2;^)/3(Aa::^)^ — 1, which 
measures how the shape of the distribution of Aa; differs 



from a Gaussian, decreases toward zero as At increases. 
(The data fluctuates somewhat for = 1.38cZ/sec, pre- 
sumably due to intermittency.) This suggests a transition 
from super to normal diffusion. 

As shown in Fig. |3Ia)-(b), the scaling of the mean- 
square displacements does, in fact, change from super- 
diffusive, (Aa;2) cx A^i-^ and (Az^) oc At^-^ to diffu- 
sive, (Aa;^) oc (Az^) oc Ai. The normal diffusion in 
long time scales is consistent with previous studies of 
dense drainage where particles were tracked with lower 
time resolution [i^, I^. Curiously, the super-diffusion 
is slower than in the usual case of ballistic transport 
in fluids, (Ax^) oc (Az^) oc Ai^, where particles move 
freely between collisions. Sub-ballistic scaling and non- 
Gaussian statistics at short times suggest that dense 
granular flows differ from classical fluid flows. 

The differences become even more obvious upon vary- 
ing the flow rate. In a fluid, this causes a linear in- 
crease in Pe because the mean flow has no affect on 
molecular diffusion due to thermal fluctuations. Here, 
as shown in Fig. |3Ic)-(d), the measured diffusion co- 
efficients, e.g. Dx, = limAt^oo(Ax^)/2At, are actually 
proportional to the flow speed (with Dz ~ 2.1Dx, consis- 
tent with the discussion above), so the Peclet numbers, 
Pe^ = Vzd/Dx ~ 321 and Pcz — Vzd/Dx ~ 150, are 
roughly constant. This suggests that diffusion and ad- 
vection are caused by the same physical mechanism, such 
as a passing void. The measured Peclet numbers, how- 
ever, are two orders of magnitude larger than predicted 
by the void model. 

Since D^-, oc v^, we plot the mean square displace- 
ments versus the mean distance dropped in the labora- 
tory frame, v^At. Remarkably, as shown in Fig.|2ta)-(b), 
this collapses all of our data for different flow speeds onto 
a single curve, not only in the diffusive regime, but also 
in the super-diffusive regime. (The data for the smallest 
flow speed again differs somewhat.) A smooth crossover 
from super to normal diffusion occurs after particles have 
fallen roughly one particle diameter. 

Although advection dominates particle dynamics 
(Pe » 1), diffusion causes a gradual rearrangement of 
the 'cage' of nearest neighbors. To investigate this mix- 
ing directly, we measure the topological correlation func- 
tion, C(At), defined as the fraction of nearest neighbor 
pairs preserved from times t to t + At, averaged over all 
t. We chose the cutoff for a nearest neighbor, l.Sd, as the 
first minimum of the radial distribution function, (which 
yields coordinations near 0.59). As shown in Fig. 0Ja), 
the data for C(At) collapses when plotted versus v^ At, in 
the sense that no systematic dependence on is observed 
(except perhaps for the smallest orifice widths), so in 
Fig. EJb) we plot the average over all experiments in the 
continuous flow regime {W > 16 mm or v^ > 5.59 d/ sec). 

The cage correlation function in Fig. Efb) exhibits 
a clear crossover, which closely parallels the ones for 
mean displacements in Fig. |3Ia)-(b). In the superdif- 
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FIG. 4: (a) The topological cage correlation function, C{At), 
versus the mean distance fallen, VzAf, and (b) the average 
over all experiments in the continuous-flow regime {W > 
16 mm), compared with an exponential fit in the diffusive 
regime (dotted line). 



fusive regime, C(At) decreases fairly quickly (with a 
decay length of roughly 20d), but after falling more 
than one particle diameter the rate of decrease (neigh- 
bor loss) slows considerably. Since the topology re- 
mains more than 90% intact within the observation 
window, the precise form of the long-distance decay is 
uncertain, but a least-squares exponential fit, (C) ~ 
0.976 exp(—V2Ai/200c?), yields a "cage breaking length" 
of 200d. This gives direct evidence that the flow is char- 
acterized by long-lasting contacts, as is obvious to the 
naked eye. It also firmly rejects the void model because 
any void-particle exchange removes roughly half of the 
neighbors of a particle as it falls by only one diameter. 

To counter the argument that a coUisional regime may 
exist below the experimental resolution {At ^ 1 ms), 
we show that this is inconsistent with the fact that diffu- 
sion and mixing depend only on geometry (Figs. 1213 . In 
the standard model of a coUisional g particle drop- 
ping a distance, L, experiences an average of N collisions 
which must dissipate its gravitational potential energy: 
m g L = (1/2) (1 — e^) to v^, where m is the mass, g the 
gravitational acceleration, e the restitution coefficient, 
and Yr the mean relative velocity. To be consistent with 
our data, N should depend on L, but not the flow speed, 
Vz- Although V,. is unknown, we can make two estimates 
— both of which lead to a contradiction. The first starts 
from the natural formula, v^ « {{Ax^) + (Az^)) / At'^ with 
fixed At = 1 ms, which suggests v^ oc v^-® by looking at 
the initial slope of (Az^) in Fig.|3|^b). The second follows 



would typically correlate with v^, so velocity-dependent 
restitution cannot compensate for the changes in v^.) 

More generally, in slow granular flows it seems that 
"granular temperature" may not be a useful concept. 
Figures |3| and 0] clearly show that fluctuations depend 
only on the distance fallen, and yet any notion of tem- 
perature should increase with the flow speed. The fact 
that the nearest-neighbor topology persists for distances 
comparable to the system size also seems to cast doubt 
on the assumption of ergodicity. 

Instead, our experimental data suggests that cage re- 
arrangements are caused by the relaxation of contact net- 
works, as are believed to occur in Couette cells Such 
networks could absorb potential energy via rolling and 
sliding neighbors. The breaking of a contact could cause 
non-Gaussian fluctuations and small-scale superdiffusion 
among the particles in a network, while the gradual de- 
struction of a network (and reformation of a new one) as 
a particle falls farther than its own diameter could cause 
the observed transition to normal diffusion. All of these 
effects are dominated by the geometry of random close 
packings, which is not fully understood, even without any 
dynamics 0. 

In summary, we have experimentally investigated par- 
ticle dynamics in dense granular flows — as they occur 
in silo drainage. Consistent with the void model, we 
observe diffusion after drainage by more than a parti- 
cle diameter and Peclet numbers which are independent 
of the flow rate, suggesting that advection and diffusion 
have the same physical source. The Peclet numbers and 
cage-breaking lengths, however, are much larger than 
predicted, so the question of an appropriate microscopic 
model is left open. 
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of Energy (grant #DE-FG02-02ER25530), the Norbert 
Weiner Research Fund, and NEC Fund at MIT and by 
the National Science Foundation (grant #DMR-9983659) 
at Clark. 
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from direct measurements |22l | of v^ yielding v,. cx y 
either case, N would not be constant. (Note that (1- 



In 

,2\ 



[1] 

[2] 

[3] 
[4] 
[5] 
[6] 

[7] 

[8] 

[9] 



H. M. Jaeger, S. R. Naegel, and R. P. Behringer, Rev. 
Mod. Phys. 68, 4, 1259 (1996). 

J. M. Ottino and D. V. Khakhar, Annu. Rev. Fluid Mech. 
32, 55 (2000). 

D. Talbot, MIT Technology Review 105, 54 (2002). 
J. Litwiniszyn, Bull. Acad. Pol. Sci. 11, 593 (1963). 
J. MuUins, J. Appl. Phys. 43, 665 (1972). 
R. M. Nedderman and U. Tuzun, Powder Tech. 22, 243 
(1979); 23, 257 (1979). 

J. R. Prakash and K. K. Rao, J. Fluid Mech. 225, 21 
(1991). 

R. M. Nedderman, Statics and Kinematics of Granular 
Materials, 241 (Cambridge University Press, 1992). 
J. Mullins, Powder Tech. 9, 29 (1974); 23, 115 (1979). 



5 



[10] A. Samadani, A. Pradhan, and A. Kudrolli, Phys. Rev. 

E 60, 7203 (1999). 
[11] H. Caram and D. C. Hong, Phys. Rev. Lett. 67, 828 

(1991). 

[12] M. Z. Bazant, http : //arxiv . org/alDs/cond-mat/0307379. 

[13] A. Medina, J. A. Cordova, E. Luna, and C. Trevino, 
Phys. Lett. A 250, 111 (1998). 

[14] S. F. Edwards, Granular Matter: An Interdisciplinary 
Approach, Ed. by A. Mehta, 121 (Springer, 1994). 

[15] A. Barrat, J. Kurchan, V. Loreto, and M. SelUtto, Phys. 
Rev. Lett. 85, 5034 (2000). 

[16] H. A. Makse and J. Kurchan, Nature 415, 614 (2002). 

[17] Z. Zhong, J. Y. Ooi, and J. M. Rotter, Engineering Struc- 
ture 23, 756 (2001). 



[18] A. ,J. Liu and S. R. Nagel, Nature 396, 21 (1998). 
[19] E. R. Weeks and D. A. Weitz, Phys. Rev. Lett. 89, 
095704 (2002). 

[20] S. S. Hsiau and M. L. Hunt, J. Fluid Mech. 251, 299 
(1993). 

[21] V. V. R. Natarajan, M. L. Hunt, and E. D. Taylor, J. 

Fluid Mech. 304, 1 (1995). 
[22] N. Menon and D. J. Durian, Science 275, 1920 (1997). 
[23] D. Howell , R. P. Behringer, and C. Veje, Phys. Rev. 

Lett. 82, 5241 (1999). 
[24] S. Torquato, T. M. Truskett, and P. G. Debenedetti, 

Phys. Rev. Lett. 84, 2064 (2000). 



